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Here we summarize our recent results of high-resolution computer simulations on the turbulent 
amplification of weak magnetic seed fields showing that such fields will be exponentially am- 
plified also during the gravitational collapse reminiscent to the situation during primordial star 
formation. The exponential magnetic field amplification is driven by the turbulent small-scale 
dynamo that can be only observed in computer simulations if the turbulent motions in the cen- 
tral core ai'e sufficently resolved. We find that the Jeans length, which determines the central 
core region, has to be resolved by at least 30 grid cells to capture the dynamo activity. We con- 
clude from our studies that strong magnetic fields will be unavoidably created already during 
the formation of the first stars in the Universe, potentially influencing their evolution and mass 
distribution. 



1 Introduction 

Magnetic fields are ubiquitous in the local Univers^and there is growing evidence of their 
presence also at high redshift^EEl The seeds for these fields could be a relic from the 
early Universe, possibly arising due to inflation or some other phase transition procesJ^. 
Alternatively, they could be generated by the Biermann battery or the Weibel instabil- 
itjUE. Regardless of their physical origin, most models predict weak field strengths and/or 
have large uncertainties '° Consequently, magnetic fields are thought to be irrelevant for 
the formation of the first stars and galaxies—; Here we show that strong and dynamically 
important magnetic fields will be generated from weak initial magnetic seed fields in the 
presence of the small-scale turbulent dynamo during the collapse of primordial halo£3Ell 



The properties of the small-scale dynamo have been explored both in comput er simula- 
tions of driven turbulence without self-gravity and in analytic model£iniG3ill as well as 
in the context of magnetic fields in galaxy clusters'^. Analytic estimates show that the 
small-s cale dyn amo could be important already during the formation of the first stars and 
galaxieslEOH 

Field amplification via the small-scale dynamo requires both turbulent gas motions and 
high magnetic Reynolds numbers. The presence of such turbulence is suggested by cos- 
mological hydrodynamical simulations of first star formation^^ where it also plays an 
important role in regulating the transport of angular momentum. High magnetic Reynolds 
numbers are expected in primordial gas as it follows closely the conditions of ideal mag- 
netohydrodynamics (MHD). 



2 Numerical method and initial conditions 



In this study, we focus on the gravitational collapse and magnetic field amplification of the 
inner parts of a contracting primordial halo. The initial condit ions for our computer sim- 
ulation were motivated from larger-scale cosmological model^mmH. Thus, we set up a 
super-critical Bonnor-Ebert (BE) sphere with a core density of pbe — 4.68 x 10~^° g cm~^ 
('^BE — 10^ cm^'^) and a small amount of rotation, £'rot/-E'g = 0.04. We use a random 
initial velocity field with transonic velocity dispersion, and a weak random magnetic field 
with Srms ~ InG. Both the turbulent energy and magnetic field spectra were initialized 
with the same power law dependence, cx fc^'^, with most power on large scales. Consis- 
tent with previous studies that follow the thermodynamics during the collapse, we adopt 
an effective equation of state with 7 = dlogT/dlogp +1 = 1.1 for all densities relevant 
to this studj^*" For this setup, we solve the equations of ideal magnetohydrodynamics 
(MHD) including self-gravity with an adaptive-mesh refinement technique^**, which guar- 
antees that the critical length scale of gravitational collapse is always resolved with the 
same number of cells regardless of density. This scale, the local Jeans length 

1/2 



where Cg and G are sound speed and gravitational constant, is set by the competition be- 
tween gravity and thermal pressure. This is also the turbulent injection scale below which 
the small-scale dynamo is active. We note that the efficiency of the dynamo process de- 
pends on the Reynolds number and is thus related to how well the turbulent motions are 
resolvecffS. Higher resolution results in larger field amplification. To demonstrate this ef- 
fect, we perform five computer simulations where we resolve Aj by either 8, 16, 32, 64 or 
128 cells. Our results indicate that a minimum resolution of 30 cells per Jeans length is 
required to see an exponential growth of the magnetic field. 



3 Results 

3.1 Physical properties 

We find that the dynamical evolution of the system is characterized by two distinct phases. 
First, as the initial turbulent velocity field decays the system exhibits weak oscillatory 
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Figure 1. Two-dimensional slices thi'ough the center of the collapsing core at the time when the initial field 
strength has increased by a factor of ^ 10^, showing the central region of about 0.02 X 0.02 pc^ in size for our 
highest-resolution simulation (Aj resolved by 128 cells). The circle indicates the central averaging volume Vj. 
The left image shows the density and the velocity component in the xy-plane, indicating radial infall in the outer 
regions and turbulent motions in the inner core. The right image depicts the total magnetic field amplitude and 
direction. [Figure taken fronii^ 



behavior and contracts only slowly. Soon, however, the run-away collapse sets in. Figure[T| 
shows a snapshot of the central region of the collapsing core in our highest-resolution 
simulation at a time when the central density has increased by a factor of ~ 10^. The 
magnetic field strength has grown by a factor of 10^, reaching peak values of about 1 mG. 
The left image shows density and velocity structure, and the right image shows magnetic 
field strength and morphology. 

To understand the behavior of the system more quantitatively, we need to account for 
its dynamical contraction. First, we note that the physical time scale becomes progressively 
shorter during collapse. We therefore define a new time coordinate r, 

J dt/tsit), (2) 

based on the local free-fall time, tff{t) — -^/Svr/ (32G(0,„(t)), where Pmit) is the mean den- 
sity of the contracting central region. We can also define a critical volume for gravitational 
collapse, the so-called Jeans volume, Vj = 47r(Aj/2)'^/3, with Aj given by equation ([T]). 
We obtain all dynamical quantities of interest as averages within the central Jeans volume. 
This approach ensures that we always average over the relevant volume for collapse and 
field amplification. 

Gravitational compression during the collapse of a primordial gas cloud can at most 
lead to an amplification of the magnetic field strength by a factor of ^ p^^^ in the limit 
of perfect flux freezing (i.e., ideal MHD). A stronger increase implies the presence of an 
additional amplification mechanism. Starting from an initial field strength of ^ InG, our 
simulations show a total magnetic field amplification by six orders of magnitude, leading 
to a field strength of about ^ 1 mG for the case where we resolve the local Jeans length 
by 128 cells. This is illustrated in Fig. |2^. Fig.jSJ) shows that the obtained field ampli- 
fication is indeed stronger than what is expected for purely adiabatic compression, which 
demonstrates that the small-scale turbulent dynamo provides significant additional field 
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T = /dt/t„(p„(t)) 

Figure 2. Evolution of tlie dynamical quantities as a function of r = J dt/tg, defined in equation d2| for five 
runs with different number of cells to resolve the local Jeans length. Panel (a) - the rms magnetic field strength 
B rmsi amplified to 1 mG from an initial field strength of 1 nG, (b) - the evolution of -Brms/Pm showing the 
turbulent dynamo amplification by dividing out the maximum possible amplification due to pure compression of 
field lines, (c) - the evolution of the mean density pni and (d) - the rms velocity Vrma ■ The onset of runaway 
collapse commences at about r ~ 6. [Figure taken from— ] 



amplification over compression. A closer comparison of Fig. [2^ and [2J) shows that the 
amplification of the field by compression and by the turbulent dynamo are roughly com- 
parable. The time evolution of the central density is depicted in Fig. [2]:, while Fig.|2}l 
shows the corresponding rms velocity. The presence of turbulence delays the collapse until 
r ~ 6. During this time, the mean density shows some oscillations while the rms velocity 
decreases as the turbulence decays. A comparison between Fig. ^ and Fig. |2jl indicates 
that dynamo amplification takes place throughout the entire duration of the simulation, that 
is during the initial phase of turbulent decay, i.e., for r < 6, as well as during the run-away 
collapse phase (r > 6). We note that the time constant for the field amplification is roughly 
the same in both regimes. 
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Figure 3. Evolution of the dynamical quantities in the case of stronger initial fields (B = 10 ^ G) within the 
central Jeans volume as a function of t (left) and the mean density pm (right) for runs with different initial Mach 
numbers. Shown are the rms magnetic field strength Brms, the ratio Brms/Pnl^, the mean density (only in 
the left panel), the rms velocity Vrme and the ratio of magnetic to kinetic energy, Em/E]^- AH the simulations 
correspond to a resolution of 128 cells per Jeans length. [Figures taken fronC^l] 



3.2 Implications for numerical resolution 

We also point out, that sufficient numerical resolution is a crucial issue when studying 
the small-scale turbulent dynamo. We see that a minimum of 30 cells per Jeans l ength 
is required to unambiguously identify the exponential growth of the magnetic field SOU 
Seeing the dynamo process is computationally extremely demanding, and even our 128^ 
high-resolution run is by no means numerically converged. This reflects the fact that no nu- 
merical scheme to date is able to reach the enormous Reynolds numbers, which determines 
the dynamo growth rate, of star-forming, turbulent gas. Therefore, this simulations show 
no signs of saturation, thus we simply stop the calculation when the numerical cost be- 
comes prohibitively high. In our subsequent stud> we investigate the saturation behavior 
when back-reactions of the Lorentz force become important. Here we find typical values of 
£'mag/^'kin of 0.2 — 0.4 (sec Fig.js]) in agreement with calculations of non-selfgravitating 
MHD turbulence'^'"'^'. The physical dissipation scales are much smaller than the Jeans 
length and thus the growth rates obtained in our simulations are lower limits on the phys- 
ical growth rates. In Fig. |4] we visualize various quantities of our 128'' resolution run to 
give an impression of the complex structure within the central core. 
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Figure 4. (a) Spherical slice of the gas density inside the Jeans volume at r = 12 for our run with 128 cells 
per Jeans length, (b) Velocity streamlines on a linear color scale ranging from dark blue (0 km s~^) to light 
gray (5 km s~^). (c) Magnetic field lines, showing a highly tangled and twisted magnetic field structure typical 
of the small-scale dynamo; yellow: 0.5 mG, red: 1 mG. (d) Four randomly chosen, individual field lines. The 
green one, in particular, is extremely tangled close to the center of the Jeans volume, (e) Contours of the vorticity 
modulus, showing elongated, filamentary structures typically seen in subsonic turbulence, (f) Spherical slice of 
the divergence of the velocity field, white: compression, red: expansion. [Figure taken fronii^ 
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4 Conclusions 



Taken all together, our results strongly indicate that dynamically important magnetic fields 
are generated during the formation of the first stars. This has important consequences 
for our understanding of how the first stars form and how they influence subsequent cos- 
mic evolution. We know from modeling galactic star-forming clouds that the presence of 
magnetic fields can reduce the level of fragmentation, and by doing so strongly influences 
the stellar mass spectrum-*^. Furthermore, we also know that dynamos in accretion disks 
produce jets and outflow^, which remove a signicant fraction of the mass and angular 
momentum. Again, this influences the stellar mass spectrum. There are first attempts to 
study this process in the context of first star formation^"*, but more sophisticated initial con- 
ditions and more appropriate magnetic field geometries need to be considered. Once the 
first stars have formed, they are likely to produce a copious amount of ionizing photons, 
which drive huge HII regions, bubbles of ionized gas, expanding into the low-density gas 
between the halos. These dynamics could be substantially different if magnetic outflows 
drive a cavity-wave into the surrounding gas. The magnetic field may further affect fluid 
instabilities near the ionization front. 

The mechanisms discussed here are likely to work not only during the formation of the 
first stars, but in all types of gravitationally bound, turbulent objects. Highly magnetized 
gas is thus expected already in the first galaxies, which subsequently needs to become more 
coherent to form the typical field structures that are observed in the present-day Universe. 
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